Non-Gaussianity with Lagrange Multiplier Field in the Curvaton Scenario 
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In this paper, we will use SAf- formalism to calculate the primordial curvature perturbation for 
the curvaton model with a Lagrange multiplier field. We calculate the non-linearity parameters Jnl 
and gML in the sudden-decay approximation in this kind of model, and we find that one could get a 
large non-Gaussinity even if the curvaton dominates the total energy density before it decays, and 
this property will make the curvaton model much richer. We also calculate the probability density 
function of the primordial curvature perturbation in the sudden-decay approximation, as well as 
some moments of it. 
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Inflation has been remarkably successful in explaining the properties of the universe and the origin of the primordial 
perturbation which is thought of as the seed of the large scale structures. So far, there still a lot of discussions 
and works on inflation, such as 0. A single- field inflation predicts a nearly Gaussian distribution of the primordial 
power spectrum Q . On the other hand, multi- field models of inflation can lead to a large deviation from the Gaussian 
distribution, which may be observed in the future observations In fact, the multi-field models generate the non- 
Gaussianity due to the non-trivial classical dynamics on super horizon scales. Since the gravitational dynamics could 
introduce significant non-linearities that would contribute to the final non-Gaussianity in the large scale of CMB 
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I. INTRODUCTION 



anisotropics, the CMB non-Gaussianity opens a window to probe the physics of the early universe. 



There are many mechanisms to generate a large local-type non-Gaussinities, and one of them is the curvaton scenario 
In this kind of model, there would be another, weakly coupled, light inhomogeneous scalar field called curvaton, 
whose energy density could be neglected during inflation, while the early Universe is dominated by inflaton. After the 
end of inflation, the energy of inflaton converted into radiations and the Hubble parameters decreases. During the 
evolution of the curvaton, its energy density goes like oc a -3 , which increases with respect to that of radiations oc a~ 4 . 
Therefore, the curvaton can dominate the energy density of the Universe later. When the Hubble parameter becomes 
the same order of the curvation decay rate, the energy of the curvaton would be converted into radiations. Finally, 
the curvation is supposed to completely decay into thermalized radiations before primordial nucleosynthesis, in the 
meanwhile, the perturbations of the curvaton become the final adiabatic curvature perturbations that seed the matter 
and radiation density fluctuations observed in the Universe. This kind of non-Gaussianity can be described by some 
non-linearity parameters /jvx, 9nl, etc. defined below. For recent progress on the curvaton model, see Ref. [a] J?} [ij] . 
At first, we expand the curvature perturbation as 



c(t,x) = ci(t,x) + 5^-c»(t,x), (i) 



n 

n=2 



where the probability density function (pdf) of the first order term £i is Gaussian, while the higher order terms give 
rise to a non-Gasussian pdf of the full £. As usual, the non-linearity parameters Jnl and g^r, are defined by 



C - Ci + IfNLCl + ^9nl(! + O(Ci) , (2) 



or, equivlently 



5(2 25C 3 



'Electronic address: fengcj@shnu.edu.cn 
^Electronic address: kychz@shnu.edu.cn 



2 



where the numerical factors arise in order to be consistent with the Bardeen potential on large scales. The upper 
bound from the WMAP 3yr data H is \ f NL \ < 114 , the bound from WMAP 5yr data [13 is -9 < f NL < 111 at 2a 
level and the constraint from WMAP 7yr data [ll[ is Jml = 32 ± 21 at la level. The correlation functions of Fourier 
transformation of £ are used to define the primordial power spectrum, bispectrum and trispectrum as 



(C(ki)C(k a )) 
(C(k 1 )C(k 2 )C(k 3 )> 

(C(k 1 )C(k 2 )C(k 3 )C(k4)> 



(27r) 3 P c (fcx)5 3 £k 

^ 71=1 

(2^) 3 S c (k 1 ,k 2 )5 3 ^k 

^ 71 = 1 

(2^) 3 r c (k 1 ,k 2 ,k 3 ) ( 5 3 ('^k 



Thus, we have 



B c (k lf k 2 ) = b -f NL 
1 8 

TcCkx.ka.ks) = -/^ 
54 

+ 25 ffAfL 



P c (A:i)-P c (A:2) + 2perms 
P c (fci)P c (fe2)Pf(|ki - k a |) + 23perms 
V ( (ki)V ( (k 2 )Vc(k3) + 3pcrms 



(4) 
(5) 
(6) 

(7) 

(8) 



In this paper, these non-linearity parameters will be calculated in the curvaton scenario with a Lagrange multiplier 
field, which is described by the following action [l2[ 



S = Id xy/—g 



K(y,X) + \(X-V(<p) 



(9) 

where the field A is a "Lagrange multiplier" without a kinetic term and the scalar field tp could be a curvaton. Here, 



(10) 



is a standard kinetic term for the field tp, K(ip, X) is arbitrary function of ip and A, and V(ip) is an arbitrary function 
of the scalar field (p. The equations of motion for A and tp are given by 



X-V(p) = 0, 



K v - V M (Jf x VV) - W v - V M (AV^( /3 ) 



0. 



(11) 

(12) 



Here and after we denote the partial derivatives by subscripts. And the energy-momentum tensor is 

T flu = (K x + X)V ft (pV v ip-Kg flv . (13) 
dt 2 — a 2 (t)dx 2 , we consider the homogeneous Lagrange 



In the spatially flat FRW universe with the metric ds 2 
multiplier field so that Eqs. (fTT]) and ([12} reduce to 



K, 



2VK XX +K x + 2\ 



V2V V2V 

And also, the energy density and pressure are given by 



<p 2 -2V( V ) = 0, 
^VKxip - 3H(K X + A) — A = 0. 



p(\,<p) = 2V(<p)(K x + A) — A , 
p(p) = K(p,V(p)). 



(14) 
(15) 



(16) 
(17) 



In this paper, we will consider two interesting cases that studied in [l2j. The first case is K v = and V = const., 
in which case, the energy-momentum corresponds to a mixture of a cosmological constant and pressureless dust, so 
we call this model the L — Xp model. From Eqs. (fl4|) and (|15p . we get 



(P : 



/ 2~V, A + 3H(K X + A) = , 



(18) 
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with solution 



<p = V2Vt, \=^_-K x ,, (19) 



where po is an integration constant. Thus, the energy density and pressure are given by 



p = — p , p = const . (20) 

or 

The other case is K = with arbitrary V, and in this model the energy density evolves exactly the same as dust or 
matter, so we call this model D — \ip model. The equations of motion for this model are 

<p = y/2V(<p) > <f + 3H l f + V v> + \- 1 \<p = 0, (21) 
and by using the relation (p = 2V(ip), we get the solution 

A = = Pa (99) 

aV 2V(ip)a 3 ' V ; 

Therefore the energy density and pressure are given by 

P=- A , P = 0. (23) 

which is exactly the behavior of the pressureless dust or the cold dark matter. It should be noticed that the expression 
of the energy density is independent of the explicit form of the scalar potential. 

It should be noticed that in a single-field inflation, the prediction of the non-linearity parameter is related to the tilt 
of the power spectrum so if a large local- type non-Gaussianity is confirmed by the future cosmological observations 
and high level data analysis, it strongly implies that the physics of the early Universe is more complicated than the 
simple single-field slow-roll inflation. In this paper, we will use cLV-formalism 13] to calculate the primordial curvature 
perturbation for the curvaton model with a Lagrange multiplier field. And we find that, in this kind of model, one 
could get a large non-Gaussinity even if the curvaton dominates the total energy density before it decays. Our paper 
is organized as follows. In Sec. UH we calculate the non-linearity parameters for the curvaton itself and in Sec. IIII1 we 
derive in the sudden-decay approximation a non-linear equation that relates the primordial curvature perturbation £ to 
the curvaton curvature perturbation ( v . Solving this equation order by order, we obtain the non-linearity parameters 
/jv l and g;si l in the sudden-decay approximation and we also compare the results with the usual curvaton model [f| • 
In Sec. IIV1 we calculate the pdf and some moments of Q in the sudden-decay approximation. In the final section, we 
will give some conclusions and discussions. 

II. NON-LINEAR PERTURBATION WITH THE LAGRANGE MULTIPLIER FIELD 

The primordial density perturbation can be described in terms of the non-linear curvature perturbation on uniform 
density hypersurface [lj] 

C(t,x) = aiV(t,x) + i f^JL, (24) 
3 J p{t) P + P 

where N — J Hdt is the amount of local expansion, and p % , p % are the local energy and pressure respectively. Here, p 
is the homogeneous energy density in the background model, while p and p is the local density and local pressure. 

A. The L — Xcp model 

In this model, we will assume that the constant pressure p is much smaller than the energy density, so the scalar 
field behaves much like the pressureless dust, thus we have the non-linear curvature perturbation on uniform-curvaton 
density surfaces is given by [3] 

C(t,x) ¥ ,=57V(t,x)+ / p-. (25) 
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Hence, the density of the scalar field tp on spatially-flat hypersurfaces is given by 

P v \sn=o = e 3C *p v . (26) 

Generally, we can expand any field 

P = <P + y2-;Sn<P, (27) 

n=l 

where <p is the homogeneous background filed. On the other sider, the quantum fluctuations in a weakly coupled 
field could be well described by a Gaussian random field [HJ]. So for such fields, one can only keep the first order 
perturbation 8±ip and the higher order perturbations 5 n ip for n > 1 that describe non-Gassian perturbations of any 
field could be neglected. However, here we also want to estimate the effect of the non-linear quantum fluctuations in 
the curvaton field and Lagrange multiplier field at Hubble exit during inflation, we will keep to the third order of the 
fluctuations as 

ip* = (f>* + Slip* + ^S 2 p* + ^5 3 p* (28) 

A* = A* + <$iA* + -5-2K + tt^A* (29) 
2 b 

where * denotes that the quantities are evaluated at the Hubble exit during inflation. Therefore, we get the density 
fluctuation of the curvaton as 

Pip* = p v * + 2V f<JiA* + -S 2 K + g^A*^j , (30) 

where p v * = 2V(Kx + A*) — p w 2V(Kx + A*). Therefore, order by order, from Eq. (121)1) we have 

SiP v * = 2V6x\* , S 2 p v * — 2V5 2 \* , 63P1P* = 2V63X* , (31) 
and we also have e 3 ^ = then 

2V 

Cv>i = x^^iA* , (32) 



2V 



C V 2 = — 82K- 3C^i = 3(ai -!)&, (33) 



where we have defined 



2V 

C V 3 = x^M* + 9(2 - 3ar)C^ = 9(2 - 3oi + h)^ , (34) 

_ p* 8 2 K pi 8 3 \* 

ai ~27(^A~F' bl ^W^WKT' [6b > 

to estimate the effect of the non-linear quantum fluctuations of the Lagrange multiplier field, and if a\ ~ and b\ ~ 0, 
this field is almost Gaussian. Using Eqs. (|32j) - (|34[) . one can express the non-linearity parameters for the curvaton 
perturbation analogous to Eq. ^ as 

/£i = §(oi-l), 9 V NL = Y (2-301 + ftx). (36) 

Here we find f^ L = —5/2 and g^ L = 25/3 for a Gaussian A field, while curvaton field itself could be Gaussian or 
non-Gaussian, and we also have the relation 10/]^ L + 3g% L = 0. It is worth to note that in this model, the fluctuation 
of the curvaton field do not contribute to the non-linearity parameters since V is a constant. 



B. The D - \tp model 



In this model, the density fluctuation is given by 

P v * = p 9 * (1 + A-^A* + V-^-SV*) 



(37) 
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where V* = V(ip*) , p^* = 214 A» 



and 



<5A* = M* + i(5 2 A* + i(5 3 A* , (38) 
2 o 



SV» = K'<W* + ^*'<W* + ^K"(^*) 2 

+ i^3V'* + ^^i^^2^ + iv;'''( ( 5iV'*) 3 , (39) 
6 2 6 



where the prime denotes the derivative with respect to p. From Eq. (I37I) . one can see that both the fluctuation of A 
and ip fields contribute to the fluctuation of the density, so in the following, we will consider two limit case: (514 = 
and (5A* = to close and open one of them and clearly illustrate the contribution of each field to the non-linearity 
parameters. In the case of (514 = , we get the the curvature perturbation as 

= (40) 
C V 2 = 1^-3^=3(^-1)^, (41) 

C V 3 = ~ + 9(2 - 3a 2 )C^ = 9(2 - 3a 2 + b 2 )C^ , (42) 

where a 2 = A*(5 2 A»/((5iA*) 2 and 6 2 = A 2 (53A*/((5iA*) 3 . Here, one can see that these equations are the same as 
Eqs. fl3"2"l) - (|3"4")) except for the definition of a, and b{. While in the case of (5A* = 0, we get 

V 

C V 2 = — <5 2 ^* + 3 (^-^72- - 1J C^i = 3 (^a 3 + -p^- - 1J C^i , (44) 

V' / V"V V"'V 2 \ 

C^3 = ^ 3 ^+9^2-3aa-3(l-a 3 )-^r + -^3 ! -J^i 

/ v"V V"'V 2 \ 

= 9 (2 - 3a 3 + b 3 - 3(1 - a 3 )^ + -U^j & , (45) 

where 

_ % 5 2 <p>* _ K 2 5 3 y* 
° 3 — 777 Tx 72 ' ° 3 _ 7777, Tx 73 ' 

Then, using Eqs. (|43[) - (|45[) . one can express the non-linearity parameters for the curvaton perturbation analogous to 
Eq. © as 

Now, let's forget about the a 3 and 63 for a while, and take the function V as the power of p, namely, V ~ <p n , then 
the non-linearity parameters becomes 

f ^ L = ~2^' 5 ^ = 3^' (48) 

which could be larger when n is small. And here, if we take V ~ e mv , these parameters would be vanished. With 
Eq. (|48|) . one can find the following relation 

10/^ + 37^=0. (49) 

It should be noticed that, the result Eq. (|4"5)l is evidently distinguishing with that in Ref. [l6j], in which the author 
considered a simple power-law potential of the curvaton, i.e. V ~ (p n , but the equation of state of the curvaton scaclar 



depends on the value of n as w = p/p = (n — 2)/(n + 2) in their situation. However, in our model, the pressure is 
always vanished. In a special case n = 2, both our result and that in Ref. are coincident with the result for the 
quadratic potential of the curvaton Of course, one can take the other forms of the function V, like V ~ cosh(m<£>) 
to generate large non-Gassuaianity even when the curvaton density dominates the universe, as long as imp* <C 1: 



if 
■ ' . 



NL 



9nl — 



and one can also find the relation 



2sinh 2 (m<^ il ,) 2(mip*) 2 ' 
25 25 



3sinh 2 (m<^) 3(m(p*) 2 ' 



(50) 
(51) 



NL 



35 



NL 



(52) 



III. SUDDEN-DECAY APPROXIMATION AND THE PRIMORDIAL CURVATURE PERTURBATION 



We assume that the curvaton decays on a uniform-total energy density hypersurface and thus from Eq. (|24p the 
perturbed expansion on this hypersurface is SN = £, where C is the total curvature perturbations at curvaton decay 
surface. And on this surface, we have 



p r (t d , x) + p v {t d , x) = p(t d ) 



(53) 



Here, we assume that all the curvaton decay products are relativistic, then Q is conserved after the curvaton decay. 
The local curvaton and radiation density on this decay surface may be inhomogeneous and they have a conserved 
curvature perturbations when they do not have interactions with each other fl4j | : 



1 



dpi 



Q = SN+ . . 

I & Pi +Pi 

Thus, the curvature perturbation related to radiations (p r = p r /3) is 



(54) 



Cr = C - 



4 VPr 



or p r = p r e 4 ^ r . 



(55) 



and 



( v = C + 3 In 



or p r = p r e 3< - cp c) 



(56) 



Therefore, from Eq. ([53]) we have the following relation 



(57) 



where fl Vid — Pip/(p~r + P<p) is the dimensionless density parameter for the curvaton at the decay time t d . Here, we take 
the sudden decay approximation for the curvaton, namely, the curvaton particles instantaneous decay into radiations. 
And for simplicity, we will ignore the small curvature perturbation in the radiation fluid before the curvaton decays, 
i.e. Cr = 0. Then, order by order from Eq. (|57p. we have 



Ci 

C2 



fdC<pl > 
3 (1 

Td U 

54 (g NL , /. 

n 



■ ff 

: JNL 



25 



NL 



fd 



18 

fd 



r J NL 



1 8 

1 ) -y/^-4 + 10/ d + 3/J 



C 6 

Sl ) 



(58) 
(59) 

(60) 



where 
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is called the curvature perturbation transfer efficiency, see Ref. |5j. Therefore, the non-linearity parameters f^L and 
qnl are given by Eq. © as 

(62) 

5 / ^ + l( 10/d + 3/ '~ 4 )- (63) 

In the limit of fd — > 1, namely, the curvaton dominates the total energy density before it decays, we recover f^L —> ff^L 
and c/n l ~> 9nl- While, in the limit of fd — > 0, one can get large non-Gaussianity as ordinary curvaton models, but 
this is not the only way to get large non-Gassianity in the D — \ip model, because one can get large non-linearity 
parameters of the curvaton itself f NL as we mentioned before. 



9nl 



£ ( fflL 



1 

fd 
1 

7! 



9nl 



5 
2 

5f* 



5(2 + f d ) 



NL 



25 

IT 



25 
Wd 



rip 

J NL 



A. The L — Xip model 



By using Eqs. flU, ([62]) and (|63]), we get 



fNL — 777 



9NL 



5ai 
256i 



5(2 + f d ) 
6 

25ai 



25 ai , 25(10/ rf + 3/ d 2 ) 



6/J 3/ rf 6 



54 



125 
~5T 



(64) 
(65) 



for the L — Xip model. Then, if fd <C a\, f^L could be large, but g^L could be large or small depending on the exactly 
values of Oi, b\ and fd- Otherwise, if ai, b\ are much smaller than fd, the terms with the factor a\ or b\ could be 
neglected, then 



NL 



(. 



fd , 9nl = 25 



5fd 
27 



fj 
18 



o 

54 



(66) 



Thus, we have —5/2 < f^L < —5/3 and 125/54 < g^L < 25/3 for a Gaussian A held. And also we have the relation 
25/tvl + l&gNL — in the limit of fd — > in this situation. 

To illustrate the above analytic study clearly, we plot the non-linearity parameters fNL and g^L as the function 
of f d (the transfer efficiency) in Fig. [T] and Fig. [5J and we also plot the results in a usual curvaton model [|| without 
nonlinear evolution of the curvaton field between the Hubble exit and the moment of its decay with a black dashed 
curve in these figures. 
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FIG. 1: The non-linearity parameters /jvl (left) and gNL (right) as a function of the curvature perturbation transfer efficiency 
fd — Ci/Cpi i Q the L — \ip model and the usual curvaton model (black dashed curve). The red solid curve corresponds to 
ai = 0.4, while the blue dot-dashed curve corresponds to ai = 0.6 and here we have set bi = 0. If ai = 0.5, &i = 0, then the 
curves predicted in the L — Xip model and the usual curvaton model are the same. 



From Fig. Q] , one can see that the shape of the curves in the L — Xip model are almost the same as that in the usual 
curvaton model. They would coincide together when a\ = 0.5, b\ = 0. If b\ does not vanish, see Fig. [21 the shape of 
the curves of gjq j, in the L — Xip model are much more different with that in the usual curvaton model when fd^O 
due to the contribution of the first (~ f^ 2 ) terms in Eq. (|6"5")) . 
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FIG. 2: The non-linearity parameters g^L as a function of the curvature perturbation transfer efficiency fd = £1/(^,1 with 
a — 0.4(left) and a — 0.6(right) in the L — \ip model and the usual curvaton model (black dashed curve). The red solid curve 
corresponds to 61 = 0.05, while the blue dot-dashed curve corresponds to bi =0.1 in both figures. 



The D - A</j model 



In the case of 5V* = 0, the result is the same as that in the L — Xip model except for the definition of dj and bi 
we will focus on the case of <5A» = in this model. By using Eqs. (|4"7) . ((521) and (1551) . we get 



Inl 



9nl 



5 

W 

25 

w d 



f'3 



vi'v* 

da 



K 2 



5(2 + f d ) 
6 

I * * 

~vp~ 



25 
3^ 



25 
~6~ 



a-3 



VP 



+ 



25(10/ d + 3/J) 



Now, lets forget about the 03 and 63 for a while, and take the function V ~ ip" 



54 

then we have 



125 
~54~ 



Tnl 



9nl 



5 
Wd 

25 



1 



(n - 2)(n - 
6fjn 2 



:fd , 

1 

3^ 



5id 
27 



18 



5 

54 



so 

(67) 
(68) 

(69) 
(70) 



which could be larger when n is small. Of course, one can also get large non-Gaussianity with fd <C 1 when 1. If 
n = 1, then /jvx and <7./vl have the same value as that in the L — \tp model, see Eq. (I66[) . If n = 2 and in the limit 
of fd — > 0, we have f^L — > 5/(4/^), g^vL — > — 25/(6 fd) and the relation 10 fNL + ^9nl — 0, which means gjvL has the 
same order as /jvl- However, if n 7^ 1, 2, we have 



fNL 777- 



1 



in the limit of fd 0. And, this time 



9NL 



9nl 



25(ra- 2)(ra- 1) 



6/Jn 2 



2(n 



3(n - 1) 



2) f2 



7VL 5 



(71) 



(72) 



which is much larger. The non-linearity parameters /jvl and g^L as the function of fd are plotted in Fig. [3] with 
V ~ <p n , and the results in a usual curvaton model with linear evolution of the curvaton field between the Hubble 
exit and the moment of its decay are plotted with black dashed curves in these figures. From Fig. [3J one can see that 
the shape of the curves are very similar except for a small value of n with red solid curves, which means one can get 
large non-Gaussianity (|/jvi| 3> 1) even if the curvaton dominates the total energy density before it decays (fd —> 1) 
by setting a small n, e.g. n = 0.1. 



If we take the function V 



Tnl 



9nl 



cosh(m<p), then we get 
5 5(2 + f d ) 



= 25 



2fd tanh (mif* 
1 



6/, 



2 tan< 



(rrup* 



6 

1 

Wd 



tanh (mp*) 



54(10/, 



3/ d 2 



(73) 
(74) 
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fd fd 

FIG. 3: The non-linearity parameters /jvl (left) and qnl (right) as a function of the curvature perturbation transfer efficiency 
fd = Ci/Cvi m the case of SX, — in the D — Ay model with V ~ y n and the usual curvaton model (black dashed curve). The 
red solid curve corresponds to n — 0.1,0.4 from top to bottom, while the blue dot-dashed curve corresponds to n — 1.5 and 
the orange dotted curve corresponds to n = 5. If n — 2, then the curves in the D — Ay model will coincide with that in the 
usual curvaton model. Here we have set a$ — 63 = for simplicity. 



As we mentioned before, when my* <C 1, /jvl and <?jvl could be larger even when fd 
we get 



/jvl 



2/^tanh (my*) 



, 9nl -> 



25 



6/J tanh 2 (my*) 



1. Taking the limit of fd — > 0, 



(75) 



and the relation 



<?ATZ = -tanh 2 (my*)/^ L . 



(76) 



The non-linearity parameters /jvl and <7jvl as the function of fd are plotted in Fig. [4] with V ~ cosh(rny). Again, 
the black dashed curves in these figures represent the results in a usual curvaton model with linear evolution of the 
curvaton field between the Hubble exit and the moment of its decay, and the shape of the curves are not similar. For 
a small value of my* = 0.1 with red solid curves in the left figure, one can get large non-Gaussianity (|/nx| S> 1) 
even if the curvaton dominates the total energy density before it decays (fd — > 1). 
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FIG. 4: The non-linearity parameters /jvl (left) and qnl (right) as a function of the curvature perturbation transfer efficiency 
fd = Ci/Cpi m t ne case °f — in the D — \ip model with V ~ cos(my) and the usual curvaton model (black dashed curve). 
The red solid curve corresponds to my* = 0.1, while the blue dot-dashed curve corresponds to my, = 0.5 and the orange 
dotted curve corresponds to rrup, — 2.0. Here we have set as = 63 = for simplicity. 



IV. PROBABILITY DENSITY FUNCTION 

In this section, we will follow the method in [n| to calculate the probability density function (pdf) of curvature 
perturbation. At first, we shall briefly review this method. Let's assume there are two random variables y and z, and 
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the functional dependence of z on y is z = z(y), which is a bijection. If the pdf of y is /(y), then the probability of z 
being in the interval {z\, z 2 ) is given by 



P(zi < z < z 2 ) 



r z 2 


dy_ 




dz 



Kv)dz . 



where the absolute value is need when y{z) is a decreasing function. Hence, the pdf of z is 

dy(z) 



/(*) = 



dz 



(77) 



(78) 



where the derivative could be replaced by the Jacobian determinant in the multi-variable case. 

Since the first order perturbation d only depends linearly on the initial Gaussian field perturbation, one can take 
Ci as a Gaussian "reference" variable with mean /i^ =0. In the sudden decay approximation, we have found an 
analytic functional dependence Q = C(Ci)j but the mapping is not always a bijection. Calling these values Cii, one can 
calculate the pdf of the non-linear primordial curvature perturbation 



/(o = E 



dCi 



fg((li) 



Ci=Cii 



where /(Ci) is the Gaussian pdf with mean \i = and variance a 1 = cr^: 

1 



/ s (Ci) 



-d 2 /(2<) 



2 ™l 



(79) 



(80) 



For simplicity, in the rest of this section, we will neglect all the non-linear fluctuation of the initial field, namely, we 
will set ai = bi = in all models. 

Actually, the non-Gaussianity could be described quantitatively by calculating the moments of the pdf and 



m z (i) = / (z - /J,) 1 j '[z)dz , 



(81) 



is called the i moment. Here, the mean can be calculated as 

M = y zf{z)dz, (82) 

where the pdf /(z) satisfies J f(z)dz = 1. Conventionally, the second moment is the variance (crj), the third moment 
is called skewness, and the fourth moment kurtosis. For a Gaussian pdf, any odd moment (with i > 3) is zero, since 
the probability density is symmetric around the mean, while the even moments could be easily calculated by partial 
integrating, e.g. m(4) = 3ct 4 , m(6) = 15cr 6 , m(8) = 105ct 8 , m(10) = 945ct 10 , m(12) = 10395ct 12 , etc., see Ref. [H. 
Any departure from these values indicates the pdf is non-Gaussian, namely, there is an asymmetric deviation from 
Gaussianity, if odd moments are not zero, and if the pdf is more (or less) sharply peaked than the Gaussian if even 
moments are smaller (or larger) than that in the Gaussian case. Therefore, the set of moments encodes the same 
information of non-Gaussianity as the fully non-linear £ or its expansion. Next, we will calculate the pdf and illustrate 
some moments of the pdf in our models. 



A. The L — \ip model 



From Eqs. j26]), J3DJ), ([32]) and ([58]), we have 



While, from Eq. ([57]) with £ r = 0, we have 



e 3C. = ! 



3Ci 
fd 



e K v = 3 + / rf c 3C + 3 f d 3 C -C 



4/d 



4/d 



(83) 



(84) 
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Thus, combing Eqs. (JS3J) and ([8l)l. we get 



and 



_1 4- 3 + /rf r 3C i 3 ^~ 3 r -c\ 

4/„ 4/ d J ' 



4 4 



Hence, the non- Gaussian probability density function for £ is 



/(0 = 



1 



2 ™i 



3 + /, 



exp 



f 1 



3 + /d , 3/d — 3 



4/d 



4/ d 



(85) 



(86) 



(87) 



15 000 



10000 



5000 




1.0010 
1.0005 

S5 i.oooo 

0.9995 



-0.0001 -0.00005 0.0000 0.00005 0.0001 

i 



0.9990 » ' 1 

-0.0001 -0.00005 0.0000 0.00005 0.0001 



FIG. 5: The L — Aip model. Left: Pdfs at fd = 0.8, (/jvl = —7/3). The black dashed curve is the pdf /, while the red solid 
curve is the Gaussian reference, f g . Right: The ratio of non-Gaussian pdf to the Gaussian one. 
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FIG. 6: Third, fourth, fifth and sixth (from left top to right bottom) moments of the pdf of the primordial curvature perturbation 
£ as a function of the linear transfer parameter, fd in the L — Xtp model. 



In Fig. 03 we compare the fully non-linear pdf f{Q to the Gaussian / s (Ci), and it shows that / is virtually 
indistinguishable from the Gaussian /„. But, we also plot f/f g which reveals the non-Gaussianity in this model. In 
Fig. [51 we plot the moments from the third up to the sixth one as a function of fd- 
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B. The D-Xip model 

In this model, we will focus on the case of 5A* = 0. From Eq. (|26|) , (|37|) and (l43l) . we have 



(88) 



and by combing Eq. (|84j) . we can get the equation for d for a given function V . Next, we will consider two types of 
V and get the non-Gaussian probability density function for £ by solving the equation as in the L — \tp model. 



1. V~<p n 



In this case we get the following equation for £, 



(89) 



where 



, i 3 + fd 3C 3/ d - 3 _ c 
a = 1 — — e s — — e ^ . 



4/ d 



4/ d 



(90) 



Here we have used Eqs. (|84p. (|88p and (j58")l . If n = 1, the results would be the same as that in the L — Xip model, 
while if n = 2, we recover the results of the usual curvaton model. For a generic n (n ^ 1,2), the above equation 
could be rewritten as 



where we have defined 



Define 



Y 3 + pY + q = , 



P 

q 



fd 3(n - 2) ' 

n 2 (n- 3) 
3(n-2) 2 (n- 1) ' 
2n 2 / d 



9(n-2)V«-l (n-l)(n-2) 3(n-2) 1 



D = 



(91) 

(92) 
(93) 
(94) 

(95) 



then, the number of real roots of the Eq. (|9"Tj) depends on the sign of D. It should be noticed that if n > 3 or n < 1, 
p is positive and then D is also positive. So, in this case, we have only one real root of the Eq. (pTT|) : 



yi=i# 8 + ul /a , u± = ^±d^ 2 



(96) 



and 



d( 



2V-D 



1/3 1/3 



(97) 



where 



Q 



! [(3 + /«i)e^ + (l-/)e-C] 



18|n- l||n - 2| 

Hence, the non-Gaussian probability density function for £ is 



(98) 



/(C) 



1 



' 1/3 1/3 I 

Q\u_ — ii_|/ 



2VD 



exp 



/d f 1/3 , 1/3 _ " 

2a 2 c T+ + "- 3(n-2) 



(99) 
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If 1 < n < 3 and n^2, there could be some regions ( £ (( a , (b), in which D < and D 



C=Ca 



D 



0. So, for 



completeness, we will give the pdf for £ in these regions in the following. If D = 0, we have three real roots of the 
Eq. (|9Tj) and two of them are equal: 



and then, we have 



(1) 



Q 



1/3 



-2/3 



£\V3 

2/ 



(2) 



(I) 



Hence, the non- Gaussian probability density function for C, is 

-2/3 

D=0 



/(C) 



27T4 



exp 



g\-2/3 
2 



exp 



2 < 



(1 



1/3 



-2/3 



1/3 



2-i 



3(n - 2) 



3(n- 2) 



And, if D < 0, we have three different real roots of the Eq. (j^T 



where 



>l „ „ /bl <?±t 

F x = 2^/ cos-, r 2l 3 = -2\/^coB— — 



Thus 



(1) 



dC, 



2 I 3 



= Q sin — sin 1 



-3/2' 



(2,3) 



Q sin sin 6M — 

^3 V 3 



Hence, the non- Gaussian probability density function for £ is 

3Q 



/(C) 



D<0 



277(7^ |psin#| 



exp 



Jl 
2al 



2\ — cos - 
V 3 3 



E 



'±7T 



exp 



3(ra-2) 

\ 2 \ 



3(n-2) 



(100) 



(101) 



(102) 

(103) 
(104) 

(105) 
(106) 



(107) 



In Fig. [71 we compare the fully non- linear pdf /(C) to the Gaussian / 9 (Ci) in the case when there is large non- 
Gaussianity. In one of them, the non-linearity parameter is very large (fd — 0.8, n = —0.001, /ml = 3126), and this 
kind of visual comparison reveals the non-Gaussianity, but in another one with (fd = 0.8, n — —0.0276, Jnl — 114), 
it shows that / is virtually indistinguishable from the Gaussian f g . But, we also plot f / f g which reveals the non- 
Gaussianity in this model. In Tab. U we give some values of the moments from the third up to the sixth one at 
f d = 0.2, 0.5, 0.8 for n = -0.001 and n = -0.0267. 



2. V ~ cosh(mi£>) 



In this case, we get the following equation for C,i from Eq. 



and 



2 tanh (m<^«) / Ci 

Jd 



2d tank 2 (mtp r ) 



= 0, 



(108) 
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FIG. 7: The D ~ Xip model with V ~ ip n . Left: Pdfs at f d = 0.8 and n = -0.001 (top, f NL = 3126), -0.0276 (bottom, 
}nl = 114). The black dashed curve ares the pdfs /, while the red solid curves are the Gaussian reference, f g . Right: The 
ratio of non-Gaussian pdf to the Gaussian one with n = —0.001 (top), —0.0276 (bottom). 



Moment 


ra= -0.0010 


n = -0.0276 


fd = 0.2 


fd = 0.5 


fd = 0.8 


fd = 0.2 


fd = 0-5 


fd = 0.8 


m(3)/cr 3 


-6.9688 x 10~ 3 


-3.5883 x 10" 3 


-2.3228 x 10 -3 


-3.5297 x 10~ 4 


-1.4023 x 10~ 4 


-8.6818 x 10~ 5 


m(4)/a 4 - 3 


-2.9970 


-2.9980 


-2.9982 


-2.9985 


-2.9985 


-2.9985 


m(5)/cr 5 


-1.2190 x 10 -3 


-5.8153 x 10~ 4 


-3.6840 x 10~ 4 


-5.5229 x 10~ 5 


-2.1891 x 10~ 5 


-1.3559 x 10~ 5 


m(6)/cr 6 - 15 


-14.999 


-14.999 


-14.999 


-14.999 


-14.999 


-14.999 



TABLE I: The value of the third, fourth, fifth and sixth moments of the pdf of the primordial curvature perturbation £ with 
n = -0.001, -0.0276 and f d = 0.2, 0.5, 0.8 in the D — \ip model with <5A, = and V ~ <p n . 



which could be rewritten as 

Z 3 +pZ + q = 0, (109) 

where we have defined 

7 Ci 1 - 2tanh 2 (m^)-l . 2 2(d-l) 2 

T A + 1' p= 3 ' 9= 27 H 9^ tanh(m^). (110) 

Thus, the solutions to Eq. f| 109[) are the same as that to Eq. (pJTl) except for the definitions of the coefficients p and 
q. Define 

~\ 3 / ~\ 2 



c=(|j+(!j. ("D 

whose sign determines the number of roots of the Eq. (|109p . It should be noticed that if tan 2 (imp*) > 1/2, p is 
positive and then D is also positive. So, in this case, we have only one real root of the Eq. (|109|) : 

ry -1/3 . ~i/3 - — 9 , ni/2 

Zi = uj_ +u_ , u± = — ± D 1 , (112) 
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and 



d(i 



2\/D 



-1/3 -1/3 
U_ — U_i_ 



where 



Q 



tanh 2 (mv5*) 
18 



(3 + f d )e^ + (1 - /)e~C 



Hence, the non-Gaussian probability density function for £ is 

1 



/(C) 



'/SI- 1 / 3 -1/3 1 
Q\u_ — M, 



2VD 



exp 



fd ( -1/3 , -1/3 1 

- 1 u + + u_ — — 



(113) 



(114) 



(115) 



If tan 2 (miy9*) < 1/2, there could be some regions of £ e (£ a , £;,), in which D < and -D|c=f a = -D|c=0> = 0- So, for 
completeness, we will give the pdf for £ in these regions in the following. If D — 0, we have three real roots of the 
Eq. (I109P and two of them are equal: 



~\ 1/3 



Z 2 =Y, 



~\ 1/3 



and then, we have 





= Q 


/ ~\ -2/3 

(?) 






_ Q 


/ -x -2/3 

(I) 


d( 






) 




2 





Hence, the non-Gaussian probability density function for £ is 



(116) 



(117) 



/(C) 



Q 



D=Q 



2 ™i 



-2/3 



exp 



-2/3 



exp 



JL 



2*1 [ 2 { 2 

2 



S\ l/ 3 



(118) 



And, for D < 0, we have three different real roots of the Eq. (|109l) : 



\p\ <? \P\ Q±k 

Zi = 2y^coB-, Z 2 , 3 = -2^cos^-, 



where 



Thus 



6 = arccos 



2 V 3 



-3/2' 



(1) 



(2,3) 



.sin^sin- 1 ^ 1 ^ 1 
3 



±7r „ : „-i2/Hp| 



Q sin — - — sin 9 — 



Hence, the non- Gaussian probability density function for £ is 

3Q 



/(C) 



D<0 



E 



2^(7^ |j5sin 
8±n 



sin - 
3 



exp 



/J LM 



2a 2 r 

Ci 



2j T COS 3-3 



exp 



2aj 

Ci 



„ \p\ 9±n 1 
2 Vf C ° S ^ + 3 



(119) 
(120) 

(121) 
(122) 



(123) 
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FIG. 8: The D — Xtp model with V ~ cosh(mvs). Left: Pdfs at f d = 0.8 and rmp, = 0.003 (top, f NL = 3472), 0.165 (bottom, 
}nl = 114). The black dashed curve ares the pdfs /, while the red solid curves are the Gaussian reference, f g . Right: The 
ratio of non-Gaussian pdf to the Gaussian one with m<p t = 0.003 (top), 0.165 (bottom). 



Moment 


mtp, — 0.03 


rrvp, — 0.165 


fd = 0.2 


fd = 0.5 


fd = 0.8 


fd = 0.2 


fd = 0.5 


fd = 0.8 


m(3)/cr 3 


1.4076 x 10~ 2 


4.3669 x 10~ 3 


2.6789 x 10~ 3 


3.5764 x 10~ 4 


1.4406 x 10~ 4 


9.0822 x 10~ 5 


m(4)/cr 4 -3 


-2.9931 


-2.9977 


-2.9982 


-2.9985 


-2.9985 


-2.9985 


m(5)/<r 5 


3.4784 x 10 -3 


7.4967 x 10~ 4 


4.3509 x 10~ 4 


5.5869 x 10" 5 


2.2490 x 10~ 5 


1.4178 x 10~ 5 


m(6)/cr 6 - 15 


-14.999 


-14.999 


-14.999 


-14.999 


-14.999 


-14.999 



TABLE II: The value of the third, fourth, fifth and sixth moments of the pdf of the primordial curvature perturbation £ with 
rrup, — 0.003, 0.165 and fd = 0.2, 0.5, 0.8 in the D — \(p model with <5A» = and V ~ cosh(m^). 

In Fig. [5J we compare the fully non-linear pdf /(£) to the Gaussian / 9 (Ci) m the case when there is large non- 
Gaussianity. In one of them, the non-linearity parameter is very large {fd = 0.8, rrup* = 0.03, Jnl = 3472), and this 
kind of visual comparison reveals the non-Gaussianity, but in another one with (fd — 0.8, nrap^ — 0.165, Jnl = 114), 
it shows that / is virtually indistinguishable from the Gaussian f g . But, we also plot f/f g which reveals the non- 
Gaussianity in this model. In Tab. UTJ we give some values of the moments from the third up to the sixth one at 
f d = 0.2, 0.5, 0.8 for rrup* = 0.03 and m<p* = 0.165 . 

V. CONCLUSION 

We have used i5A/"-formalism to calculate the primordial curvature perturbation for the curvaton model with a 
Lagrange multiplier field in two interesting cases. We have calculate the non-linearity parameters fNL and gNL in 
the probability density function of the primordial curvature perturbation in the sudden-decay approximation, as well 
as some moments of it. We find that one can get a large non-Gaussianity in this kind of model even if the curvaton 
dominates the total energy density before it decays, namely fd 1, while in the usual curvaton model with quadratic 
potential, one can only get a large non-Gaussianity when the curvaton is subdominant, namely in the limit of fd — > 0, 
see Q. It should be noticed that the isocurvature perturbations are created when the curvaton fail to dominate the 
energy density while decaying. So, it will not produce large fNL comparing to that of WMAP by taking account of 
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the constraint of isocurvature perturbations [17J- So, it seems that we have escaped the constraint from isocurvature 
perturbations and give a large Jnl comparable to the result of WMAP. So, the introduction of Lagrange multiplier 
field will make the the curvaton model much richer, e.g. it also release the form of the potential. 

Furthermore, it hase been shown that with the help of a Lagrange multiplier field, one can proposed a way to unify 
dark matter and dark energy in a single degree of freedom, see [12j . And the Lagrange multiplier modified gravity 
may lead to cyclic behavior very easily in the cyclic cosmology, and the scenario is much more realistic in the case of 
scalar cosmology [l8l ]. So, we conclude that the Lagrange multiplier field could play a very interesting and important 
role in the construction of cosmological models, and it also interesting to study the primordial nonlinear structures 
and black holes [lSl in such kind of curvaton scenario. 
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